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We present the results of high resolution absorption-line spectroscopy of 3 
face-on galaxies, NGC 98, NGC 600, and NGC 1703 with the aim of searching 
^ ! for box/peanut (B/P)-shaped bulges. These observations test and confirm, for 

the first, time the prediction that face-on B/P-shaped bulges can be recognized by 
two minima in the profile along the bar's major axis of the fourth Gauss-Hermite 
'Eo ! moment of the line-of-sight velocity distribution (LOSVD). In NGC 98, a clear 

double minimum in is present along the major axis of the bar and before the 
i end of the bar, as predicted. In contrast, in NGC 600, which is also a barred 

^ I galaxy but lacks a substantial bulge, we do not find any significant kinematic 

00 '. signature for a B/P-shaped bulge. In NGC 1703, which is an unbarred control 
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6 
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galaxy, we found no evidence of a B/P bulge. We also show directly that the 
LOSVD is broader at the location of the minimum in NGC 98 than elsewhere. 
O I This more direct method avoids possible artifacts associated with the degeneracy 

Q I between the measurement of line-of-sight velocity dispersion and /i4. 
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Subject headings: galaxies: bulges — galaxies: evolution — galaxies: formation 
— galaxies: kinematics and dynamics — galaxies: spiral 



Introduction 



Roughly one-quarter of the visi ble light emitted by s t ars in the local unive rse comes 
out of the bulges of disk galaxies (jPersic &: Saluccil Il992l : iFukugita et al.l Il998l ) . Under- 
standing how bulges form is therefore of great importance to developing a complete pic- 
ture of galaxy formation. The processes by which bulges form are still debated. On the 
one hand, the merger of dwarf-sized galactic sub units has been suggested as the main 
path for bulge formation (IKauffmann et al.l Il993l ). which is support ed by the relatively 
homogeneous bulge stellar populations of the Milky Way and M31 (IFerreras et al.l l2003l : 



Zoccali et al. 2003 



Stephens et al.ll2003l ). Bulges formed in such mergers are termed 'clas- 



driven gas inflows, bending ins 


tabilities, and clump insta 


bilities ( 


Combes & SandersI 1981: 


Pfennieer 1984: Combes et al. 


199G: Pfennieer & Friedh 


199] 


L: Raha et al. 1991: 


Noeuchi 


1999: Immeh et al. 2004: Athanassoulal 2005: 


Debattista et al. 


200( 


3). Moreover, the growth 



of the bulge out of disk mat erial may also be externally triggered by sa tellite accretion 
during minor merging events (lAguerri et al.l 1200 ll : lEliche-Moral et al.l |2006| ) or also via ex- 
ternal 'secular' proce sses. Evidence for secular bulge formation includes the near-exponential 



1996 



de Jong 



1996 



MacArthur et al 



MacArthur et al 



bulge light profiles (lAndredakis fc Sanders 



CaroUo et al. 



1998 



1994 



Andredakis et al. 



Prieto et al 



2001 



1995 



Carollo et al 



Courtea u et al.l 



2001 



20031 ) , a correlation between bulge and disk scale le ngths (Ide Jong 



Carollo 1999 



1996 



bulges and inner di sks (ITerndrup et al 



20031 : Aguerri et al. 



2005 



Mendez-Abreu et al 



1994 



Peletier fc Balcelli 



20081 ) , the similar co 



19961: 



Courteau et al. 



ors of 



1996 



20021 



Carollo et al.l l2007l ) . substantial bulge rotation (iKormendyl Il993l : iKormendv et al. 
and the presence of B/P-shaped bulges in ~ 45% of edge-on disk galaxies (jLiitticke et al. 
20001). A review of secular 'pse udo-bulge' formation and evidence for it can be found 



m 



Kormendy &: KennicuttI (120041 ). Standard c old dark matter cosmology predicts that 



galaxies without cl assical bulges should be rare (ID'Onghia fc Burkertll2004l ). Not only are 
they not rare (e.g., iLaurikainen et al.ll2007l ). they are more prevalent among normal galax- 
i es than classical bulge s, exacerbating the disagreement between theory and observations 
(jPebattista et al.l l2006l ). It is therefore important to determine which bulges are of the 
classical versus pseudo variety, and which are a mix of both. 

A/"-body simulations show that barred galaxies have a tendency to develop B/P bulges 



fICombes fc SandersI Il98ll : ICombes et all Il990l : IPfenniger fc Friedlil Il99ll : iRaha et all Il991 
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Debattista et al.l 120061 ) . The fraction of edge-on gal axies with B/P bulg e s is comparable to 



the fr action of disk galaxies containing bars (~ 60% lEskridge et al.ll2000l : iMarinova &: Jogee 



bulges has been esta 


jlished by the kinematics o 


" gas and stars ( 


Kuiiken & Merrifield"l995': 


Merrifield & Kuiiken 


1999: 


Bureau & Freeman 


1999: 


Bureau & Athanassoula 


1999. 12005.: 


Vega Beltran et al. 


1997: 


Athanassoula & Bureau 


I999I: 


Chune & Bureau 


2004 


). However, 



the degeneracy inherent in deprojecting edge-on galaxies makes it difficult to study other 
properties of the host galaxy. For example, simulations show tha t a bar can produce a B/P 
shape even if a massive cla ssical bulge formed before the disk (lAthanassoula fc Misiriotis 
2OO2I : iDebattista et al.ll2005l ). Understanding the relative importance of classical and pseudo- 
bulges requires an attempt at a cleaner separation of bulges, bars and pea nuts, which is easi- 
est to accomplish in less inclined systems. A t least in two cases, NGC 4442 (IBettoni &: Galletta 
I994J ) and NGC 7582 (iQuillen et al.l 119971 ). a B/P bulge can be recognized at inclin ation 
i < 75°, which permits also a direct identification of the bar. IDebattista et al.l (120051 ) pro- 
posed a kinematic diagnostic of B/P bulges which properly works up to a disk inclination of 
about 30° namely a double minimum in the fourth moment, /14, of the line-of-sight velocity 
distribution (LOSVD), along the major-axis of the bar. These minima occur because at 
the location of the B/P shape, the vertical density distribution of stars becomes broader, 
which leads to a double minimum in Z4, the fourth-order Gauss- Hermite moment of the 
vertical density distribution. The kinematic moment is then found to be an excellent 
proxy for the unobservable Z4. In contrast, the increase in the vertical scale-height does not 
produce any distinct signature of a B/P bulge and the vertical velocity dispersion, a^, is too 
stro ngly dependent on the radial density variation to provide a useful B/P bulge diagnos- 



tic. IDebattista et al.l (120061 ) showed that the diagnostic continues to hold even when gas is 
present since this sinks to a radius smaller than that of the B/P bulge. 

In this Letter we present the first high resolution stellar kinematics of face-on galaxies 
with the goal of search ing for B/P-shaped bulges and testing the diagnostic proposed by 
Debattista et al.l (120051 ). In ^we describe the observations of the sample and data reduction, 
while our results are presented in ^ Finally, the conclusions are given in §11 



2. Observations and data reduction 

2.1. Sample selection 

The barred galaxies NGC 98 and NGC 600 were selected in the NASA/IPAC Extragalac- 
tic Database (NED) as bright {Bx < 14) and undisturbed objects, with a disk inclination 
lower than 30° and a bar length larger than 10" to efficiently test the kinematic diagnostic. 
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a disk diameter smaller than 4' to allow an accurate sky subtraction, no strong evidence of 
dust and no bright foreground stars in the Digitized Sky Survey image. To select the sam- 
ple galaxies, their disk inclination, bar length, and bar position angle were first determined 
from ellipse fits to the 2MASS J— band archival images. The main properties of the sample 
galaxies are given in Table [H 



2.2. Long-slit spectroscopy 

The spectroscopic observations were carried out in service mode at the Very Large 
Telescope (VLT) at the European Southern Observatory (ESO) on 2005 October 2-9 and 
November 2. The Focal Reducer Low Dispersion Spectrograph 2 (F0RS2) mounted the 
volume-phased holographic grism GRIS_1028z+29 with 1028 grooves mm~^ and the O'.'T x 6'.8 
slit. The detector was a mosaic of 2 MIT/LL CCDs, each with 2048 x 2068 pixels of 15 x 15 
/im^. The wavelength range from 7681 to 9423 A was covered in the reduced spectra with 
a reciprocal dispersion of 0.858 A pixel"^ and a spatial scale of 0'.'250 pixel"^ after a 2 x 2 
pixel binning. The spectra were taken along the bar major axis of NGC 98 (P. A. = 32?2) 
and NGC 600 (P. A. = 22?1) and along the disk major axis of NGC 1703 (P. A. = 148?8). 
The total integration time for each galaxy was 3 hours, in four exposures of 45 minutes each. 
Comparison lamp exposures obtained for each observing night ensured accurate wavelength 
calibrations. Spectra of G and K giant stars served as kinematic templates. The average 
seeing FWHM was 1'.'3 for NGC 98, 0'.'9 for NGC 600, and 0'.'9 for NGC 1703 as measured 
from the ESO Differential Imaging Meteo Monitor. 

Using standard IRAF0 routines, all the spectra were bias subtracted, flat-field corrected, 
cleaned of cosmic rays, corrected for bad pixels, and wavelength calibrated as in Debattista 
et al. (2002). The accuracy of the wavelength rebinning (pa 1 kms~^) was checked by 
measuring wavelengths of the brightest night-sky emission lines. The instrumental resolution 
was 1.84 ± 0.01 A (FWHM) corresponding to o-i„st = 27 kms~^ at 8552 A. The spectra 
obtained for the same galaxy were co-added using the center of the stellar continuum as 
reference. In the resulting spectra the sky contribution was determined by interpolating 
along the outermost ~ 30" at the edges of the slit and then subtracted. 



^IRAF is distributed by NOAO, which is operated by AURA Inc., under contract with the National 
Science Foundation. 
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2.3. Photometry 



We analyzed the uncahbrated acquisition images from the VLT to derive the photometric 
properties of the sample galaxies. Isophote-fitting with ellipses, after masking foreground 
stars and bad pixels, was carried out using the IRAF task ELLIPSE. Under the assumption 
that the outer disks are circular, their inclination was determined by averaging the outer 
isophotes. All 3 galaxies have i < 30° (Tabled]). The semi-major axis length, cb, of the 
bars of NGC 9 8 and NGC 600 (Table [1]) was measured from a Fourier decomposition as in 
Aguerri et al.l (120001 ). It was calculated using the bar/interbar intensity ratio /b//ib- The 
bar intensity is defined as the sum of the even Fourier components, Iq + I2 + + 1^, while 
the inter-bar intensity is given by Iq — I2 + I4 — le,- The bar length is th e FWH M of radial 
profile of Ih/hh- This method was applied by lAthanassoula fc Misiriotij (120021 ) to analytic 
models demonstrating its accuracy in measuring ob. 



2.4. Kinematics 



The stellar kinematics of the three galaxies were measured from the galaxy absorption 
features present in the wavelength range centered on the Ca II triplet (AA 8498, 8542, 8662 A) 
using the Penalized Pixel Fitting method (pPXF. ICappellari &: Emsellemll2004l ). The spectra 
were rebinned along the dispersion direction to a logarithmic scale, and along the spatial 
direction to obtain a signal-to-noise ratio S/N > 40 A~^ It decreases to S/N = 20 A^^ at 
the outermost radii. 

A linear combination of the template stellar spectra, convolved with t he line-of-sight ve- 

l ocity distribution (LOSVD) described by the Gauss- Hermite expansion by Ivan der Marel fc Franx 
(119931 ) was fitted to each galaxy spectrum by minimization in pixel space. This allowed 
us to derive profiles of the mean velocity (fios), velocity dispersion (aios), third (h^), and 
fourth-order (h^) Gauss-Hermite moments. The uncertainties on the kinematic parameters 
were estimated by Monte Carlo simulations with photon, read out and sky noise. Extensive 
testing on simulated galaxy spectra was performed to provide an estimate of the biases of the 
pPXF method with the adopted instrumental setup and spectral sampling. The simulated 
spectra were obtained by convolving the template spectra with a LOSVD parametrized as a 
Gauss-Hermite series and measured as if they were real. No bias was found in the ranges of 
S/N and (Xios which characterize the spectra of the sample galaxies. The values of hs and /14 
measured for t he simulated spectra d iffer from the intrinsic ones only within the measured 
errors (see also lEmsellem et al.l 120041 ). 
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Results 



The photometric and kinematic profiles of the sample of galaxies are shown in Figure [TJ 
We focus here on the profiles; a full analysis of this data set will be presented in a 
forthcoming paper. For the unbarred galaxy, NGC 1703, the profile shows no sign of 
a minimum. The same is true for the almost bulgeless, barred NGC 600. In both these 
galaxies a > 30 kms~^, including within the bar in NGC 600. Thus, our failure to find an 
minimum cannot be ascribed to low spectral resolution. In NGC 98, instead, we find 
a double minimum in the /i4 profile. The tw o minima are symmetr ic with respect to the 
galaxy center (|r| ^ 5"). In the simulations of iDebattista et al.l (120051 ) the minima in /i4 are 
~ 0.05 — 0.1 deep, comparable to the minimum in the profile of NGC 98. 

We confirmed that our kinematics, and minima in h^, are not a result of details of 
the pPX F method by repeating th e analysis using the Fourier Correlation Quotient (FCQ) 
method (lBenderlll990l ) as done i n iPizzella et al.l (120081). Nonetheless expansi ons in Gauss- 
Hermite moments are degenerate (Gerhard 19931 : van der Marel fc Franx 1993) since, for the 
even moments, it is possible to obtain comparably good fits by trading uios for h^. In order 
to confirm that the signature of a B/P bulge obtained in NGC 98 is not due to a conspiracy 
of such degeneracies we introduced an LOSVD broadening measure defined for each radius 
as 

Bix,,X2) = Wix2)/Wix^), (1) 



where 



W{x) 



^max 
los 



F{v)dv, 



(2) 



mill 

los 



with F{v) the cumulative distribution of the LOSVD. The values of f™™ and v^f^ corre- 



^los '^'-'-'^ '^los 

spond to F(t>[™™) = 0.5 (1 — x) and F{v^^^) = 0.5 (1 + x), respectively. B{xi,X2) has the 
advantage of not being sensitive to the wings of the distribution, unlike the kurtosis and is 
a non-parametric measure of the LOSVD shape. Experiments with A^-body simulations, in- 
cluding those presented in IDebattista et al.l (120051 ) . showed that B{0.7, 0.9) is a useful peanut 
diagnostic in the same spirit as h^. Figured] plots B{0.7, 0.9). A clear peak in the broadening 
is apparent at \r\ ~ 5" in NGC 98. We conclude that the evidence for a B/P-shaped bulge 
in NGC 98 is robust. 



4. Conclusions 



We have identified a B/P-shaped bulge in one galaxy from a sample of three face-on 
galaxies. In the unbarred control galaxy NGC 1703 we had not expected to find a B/P 
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bulge. The failure to find a dou ble- minimum in is therefore fully consistent with previous 
results ( IChung fc Bureaij|2004l ). Of the two barred galaxies, NGC 98 has clear evidence of 
a B/P bulge while NGC 600 does not. The absence of a B/P shape in the latter galaxy is 
not surprising since it appears to not have a bulge and its bar has a lumpy structure. 



If we identify the radi us of the B/P bu 



hi, as found in simulations (Debattista et al. 



ge, -R r/p, with the location of the minimum in 
20051) ■ then we find -R r/p ~ O.SSob, where or is 



the bar semi-major axis. Similarly, iKormendy fc KennicuttI (12004 ) noted that the maximum 
radius of the boxy bulge is about one-third of the b ar radius. Simulations also produce B/P 



bulge s which are generally smaller than the bar (jShen fc SellwoodI |2004J : iDebattista et al. 
2005h . 



While observations have established the prese nce of bars in the naajori ty of systems with 
B/P-bulges, the converse has not been established. IDebattista et al.l (120051 ) demonstrate how 
the presence of a classical bulge can in some cases mask the presence of B/P bulge. Moreover, 
the ratio of radii of bulges and B/P bulges is determined by resonances. For these various 
reasons it would be very instructive to repeat measurements such as here for a sample of 
barred galaxies. 
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Table 1. Parameters of the sample galaxies. 



Object 


Type 




D25 


MO 


i 




NGC 98 


SBbc 


ISA 


1:7 


-22.4 


26° 


14" ± 1" 


NGC 600 


SBd 


12.9 


3f3 


-20.0 


21° 


12" ± 1" 


NGC 1703 


Sc 


11.9 


3:0 


-20.6 


15° 





Note. — Hubble types, apparent magnitudes, 
and diameters arc from the NASA/IPAC Extragalac- 
tic Database (NED). Absolute magnitudes are calculated 
from (NED) with distance from Vcmb (NED) and 
Hq = 75 kms~^kpc~^. Disk inclinations and bar semi- 
major axis lengths are from this paper. 




Fig. 1. — Morphology and stellar kinematics of NGC 98 (left panels), NGC 600 (central 
panels), and NGC 1703 (right panels). For each galaxy the top panel shows the VLT/F0RS2 
/—band image. The slit position and image orientation are indicated. The inset shows the 
portion of the galaxy image marked with a white box. The gray scale and isophotes were 
chosen to enhance the features observed in the central regions. The remaining panels show 
from top to bottom the radial profiles of surface brightness (extracted along the slit with 
an arbitrary zero point), i3(0.7, 0.9), /;,4, h-^, o"ios (in kms^"*^ ), and velocity v (in kms~^ , 
obtained by subtracting the systemic velocity from Vios)- The two vertical lines indicate the 
location of the minima in NGC 98. 



